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PROCEDE DE CARACTERISATION OPTIQUE DE MATERIAUX SANS 
UTILISATION DE MODELE PHYSIQUE 

DESCRIPTION 

5 DOMAINE TECHNIQUE 

La presente invention concerne un precede 
de caracterisation optique de materiaux. 

Ce procede permet de caract§riser des 
couches minces ou epaisses de ces materiaux, qui sont 
10 fornixes sur des substrats. Les grandeurs physiques, que 
ce procede permet de determiner, sont: 

- l'epaisseur d'une couche d'un mat6riau, 

- l'indice de refraction de ce materiau, et. ; 

- le coefficient d' absorption de ce yr j 
^ . $ '\ 

15 matenau. -V? 

• v 

La caracterisation optique des materiaux ,V 
est utile pour 1 1 analyse chimique de ces materiauxY 
(notamment etude des bandes d' absorption, des 
proprietes de densif ication et des proprietes 

20 d'oxydation) , dans les domaines de la raicrq- 
electronique, des capteurs, de la biologie, de la 
medecine) , ou bien pour l r analyse des epaisseurs de 
depots de ces materiaux. 

Pour des exemples d ' application, on se 

25 reportera au document [1] qui, comme les autres 
documents cites par la suite, est mentionng a la fin de 
la presente description. 

La caracterisation des proprietes optiques 
d'un matSriau est egalement utile lorsque le materiau 

30 est par la suite structure ( pour y former par exemple 
des gravures ou des rugosites) et que les proprietes de 
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diffraction optique de la structure obtenue doivent 
etre calculees (voir le document [2] ) . 

Indiquons des & present que 1" invention est 
particulierement utile lorsque la loi physique suivie 
par I'indice de refraction complexe du materiau que 
l'on veut caract^riser n'est, a priori, pas connue. 

ETAT DE LA TECHNIQUE ANTERIEURE 

On rappelle que les mesures optiques 
peuvent etre de divers es natures : 

II peut s'agir de mesures 

ref lectometriques . Dans ce cas, le coefficient de 
reflexion en intensite d l une structure est mesure sur 
un spectre (c 1 est-a-dire un intervalle) de longueurs 

d'ondes [X m , X^] . 

L 1 angle d 1 incidence de la lumiere 
d'eclairement peut etre non nul. Le coefficient de 
reflexion peut etre mesure pour plusieurs angles 
d ! incidence 9. Nous noterons R (6, X, p) le spectre 
reflectometrique, ou p est la polarisation du faisceau 
incident at X la longueur d'onde de ce dernier. 

Generalement, 1' angle 9 est nul et la 
polarisation p indeterminee . Dans le cas ou 9 n'est pas 
nul, il faut connaitre cette polarisation p. En 
general, cette derniere est de type (S) ou (P) . 

II peut s f agir aussi de mesures 
ellipsometriques . Les grandeurs me sure es sont alors les 
parties reelle et imaginaire du rapport du coefficient 
de reflexion en polarisation (P) au coefficient de 
reflexion en polarisation (S) . 
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On note generalement p= i p I exp ( j A) ce 
rapport complexe (avec j 2 =-l) . Les grandeurs 
general erne nt traitees sont Ipl, que l'on note tan(\j/) / 
et cos (A) , ou des combinaisons des deux. Par souci de 
5 generality, nous noterons S a et S 2 ces grandeurs. 

Les spectres Si, ie [1,2], sont mesurgs sur 
une plage de longueurs d'ondes t^m/ X^] . L f angle 
d' incidence peut etre quelcongue. Plusieurs spectres 
peuvent etre mesures S differents angles d 1 incidence 

10 afin d'obtenir un spectre plus riche. Nous noterons 
s (8, X) ={Si (9, X) , S 2 (6 r X)} le spectre ellipsometrique. 

De fagon complementaire, des mesures 
goniom^triques (coefficient de reflexion en fonction de 
1' angle d* incidence) peuvent etre ajoutees aux mesures 

15 • servant a la caracterisation, afin de determiner 
l'epaisseur des diverses couches, pour une ou plusieurs 
longueurs d'ondes. Ces mesures ne sont pas suffisantes 
en elles-memes puisque l'on veut determiner l'indice de 
refraction complexe sur une gamme spectrale allant de 

20 a ^m. 

Afin de simplifier l'expos6, nous noterons 
*F un ensemble de spectres ref lectometrique ou/et 
ellipsometrique (s) (et eventuellement goniometriques 
pour quelques longueurs d*ondes) . 
25 Sans perdre aucune generality, nous 

n'exposerons, dans la presente description, le mode 
d'utilisation des proced^s de.l'art anterieur et de la 
presente invention que dans le cas d'une seule couche 
mince d l un materiau, formee sur un substrat connu. 
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L'epaisseur de cette couche est notee e et 
l f indice de refraction complexe du materiau & la 
longueur d'onde X est note n* (X) . 

On rappelle & ce propos que la partie 
reelle (respectiveraent imaginaire) de cet indice de 
refraction complexe est notee n(X) ( respect ivement 
k(A,)) et appelee "indice de refraction 11 (respect ivement 
"coefficient d' extinction" ) . 

En outre on note Er (4 /<1) , *P (2) ) une fonction 
d'erreur (par exemple l'ecart quadratique moyen) entre 
deux spectres *¥ {1) et ¥ (2) . 

Par exemple, on pourra prendre, lorsque 
l'on dispose de spectres ellipsometriques sur plusieurs 
angles Q±, ie {l..n}, et d'un spectre ref lectometrique : 

+ [ji^(-i)-ji (i> a)fjflw (i) 

avec (^)={S! (1) (e i€ {l. -n}A) ,S 2 (1) (9 i€ {l. .n},X) ,R (1) (X)} 

et ¥ (2) (^-{S^ 21 (8 i€ {l. .n},X> ,B 2 m (e i6 {l. .n},X) ,R {2) (X)} 

Des facteurs de ponderation peuvent etre 

apportes £ 1 1 int^grale de fagon a ce que la fonction 

d'erreur puisse tenir compte de variations sur la 

presicion de mesure des spectres. 

La caract^risabion optique de couches de 

materiaux s'articule genera lement autour de deux 

applications : 

La premiere application est le controle 

dimensionnel du dep6t de couches minces utilisees en 

microelectronique . 
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General ement on connait bien le materiau 
depose, c'est-S-dire que l ! on connait bien 1' indice de 
refraction complexe de ce materiau aux longueurs 
d'ondes de la lumiere utilisee pour la caracterisation. 
5 Les lois suivies par 1 ' indice de refraction 

complexe sont soit tabulees, soit approch^es par des 
lois physiques connues telles que, par exemple, le 
modele de Cauchy, le modele de Sellmeier (voir le 
document [3] ) , les lois de Forouhi (voir le document 

10 [4] ) , et les lois lois d 1 oscillateurs harmoniques (voir 
le document [5] ) . Ces lois sont def inies par un nombre 
fini de paramStres . 

Par exemple, une loi de type Cauchy sans 
absorption a deux parametres est def inie de la maniere 

15 suivante : 

Re[n*(A)] =n(X)=a 0 + ^V 
lm[«*(/l)]=k(k)=0 

Lorsque 1 1 on est sur de la - valeur des 
coefficients ai(i€={l,2}) mais que 1 1 on ne connait pas 
20 l'gpaisseur, un algorithme de recherche est utilise 
afin de trouver l'£paisseur qui minimise 1 1 erreur entre 

la mesure *F et la reponse theorique ¥ compte tenu de 
I 1 indice modelis£. 

L 1 algorithme de recherche peut etre, par 
25 exemple, la methode du Simplex, la recherche Tabou, la 
methode de Levendt-Marquart ou la methode du recuit 
simul£ (voir le chapitre 10 du' document [6]). 

Lorsque 1 1 indice de refraction est 
approximatif , les coefficients ai sont integres dans la 
30 procedure d'ajustement de Y et ¥ . La recherche des 
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coefficients ai constitue un precede de caracterisation 
de 1 1 indice de refraction. 

Cependant, lorsque la loi suivie par cet 
indice de refraction est inconnue (il arrive que le 
5 materiau soit inconnu ou bien ne soit pas bien decrit 
par une loi physique connue) , ce procede reste 
approximatif et I'epaisseur risque d'etre fausse. 

La deuxieme application est la 

caracterisation des materiaux. 
10 Le procede utilise reste le meme, mis a 

part que le materiau n'est pas bien connu. C'est 
precisement la fonction d 1 indice de refraction complexe 
la plus proche de la realite qui est visee. 

Le type de loi peut etre choisi par 
15 analogie avec d'autres materiaux. Cependant, la loi 
suivie par 1 1 indice de refraction complexe peut etre 
compliquee, ce qui est par exemple le cas d'une loi 
d 1 oscillateurs harmoniques : 

[»(£) + ju E)f = i + y — ^ — 

20 Dans I 1 expression ci-dessus, j 2 =-l et 1 1 indice de 
refraction et le coefficient d 1 extinction sont exprimes 

non pas en fonction de X mais de avec E= 1240/A, 

(X en nm) . 

Dans ce cas, les coefficients des 
25 oscillateurs sont difficiles a trouver si 1 1 on n'a pas 
leur ordre de grandeur. La recherche est dif f icilement 
automat isable, les algorithraes de recherche pouvant 
donner des reponses erronees et le temps perdu peut 
etre considerable . 
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II existe une alternative a la recherche 
des coefficients : la methode point-£-point (PAP) . 
Cette methode PAP propose de ne pas choisir de loi 
physique et de rechercher 1 ! indice de refraction 
5 complexe du mater iau pour chaque longueur d T onde X±, ie 
[l...n], avec A,!=A^ et A, n =X M . 

Pour chaque X± f un algorithme de recherche 
tente de trouver l'epaisseur, 1 1 indice n(A,i) et le 
coefficient d' extinction k(Xi) qui minimisent 1 1 erreur 
10 entre la mesure *F (X±) et la r6ponse th^orique 
V a ± ,na±) ,k(^i) ,e) . 

Un tel proc^de pose un probleme parce que 
les divers points (X±, e,n(Xi) , k (X±) ) ne sont pas 
forc§ment physiquement compatibles entre eux : par'. 
15 exemple, l'epaisseur trouvee peut varier en fonction de 
la longueur d'onde et la loi suivie par 1 1 indice de 
refraction complexe, plus simplement appelee loi 
d' indice, peut presenter des discontinuites-. 

Ce procede est general ement valable 
20 seulement lorsque l'epaisseur est tres bien connue et 
que les mesures sont de trSs bonne quality . 

EXPOSE DE L' INVENTION 

La presente invention a pour but de 
remSdier aux inconv^nients precedents. 
25 Le procede objet de 1' invention permet de 

caracteriser un raateriau sans utiliser un modele 
physique, c'est-a-dire sans utiliser une loi physique 
suivie par 1' indice de refraction complexe du materiau 
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etudie. II est done tout particulierement utile 
lorsqu'une telle loi n'est pas connue . 

Ce procede constitue une alternative aux 
precedes de caracterisation connus, mentionn^s plus 
haut . II peut Stre appele "methode des nceuds" car il 
utilise des "noeuds" e'est-a-dire des points de 
coordonnees (X L n* ) , ou n* est la valeur prise par 

l 1 indice de refraction complexe a la longueur d'onde Xx 
et _i prend un nombre limit e de valeur s (entieres) . 

De fagon precise, la presente invention a 
pour objet un procede de caracterisation optique d'au 
moins une couche d'un materiau dans un intervalle de 
longueurs d'ondes [X min, X max], cette couche 6tant 
formee sur un substrat, ce procede etant caracterisS en 
ce que : 

- on effectue un ensemble de mesures de ref lectometrie 
et/ou d 1 ellipsomStrie, cet ensemble de mesures 
conduisant a un spectre mesur£, note V F, 

- on choisit m longueurs d'ondes initiales Xx »Am 
appartenant a cet intervalle, m etant un nombre entier 
au moins egal a 1, on associe, a chaque longueur 
d'onde, un indice de refraction, 

- on choisit une loi d 1 interpolation au moins pour 
1 1 indice de refraction du materiau, pour les longueurs 
d'ondes comprises entre les longueurs d'ondes 
initiales X^X^, 

- on choisit M paramdtres initiaux, M etant au moins 
<§gal a m, a savoir un indice de refraction initial n± 
pour chaque longueur d'onde initiale X±, 1 £ i £ m f 
les longueurs d'ondes initiales etant choisies de 
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maniere a pouvoir determiner par interpolation au 
moins 1 1 indice de refraction pour toute longueur 
d'onde de 1 1 intervalle [X min 7 X max], les couples (\ if 
ni) £tant appeles noeuds, 
5 - on choisit des methodes de calcul de ref lectometrie 
et d 1 ellipsom§trie, 

- on choisit aussi une fonction d'erreur Er, 
representative de l'ecart entre deux spectres *Fi et 
*P 2 , les spectres *Fi et *¥ 2 etant calcules ou mesures 

10 sur un nombre de points supSrieur au nombre m de 
noeuds, 

- k l'aide des m longueurs d'ondes initiales, des M 
parametres initiaux et de la loi d f interpolation, on 
met en ceuvre le processus d 1 optimisation suivant : 

15 — on determine un spectre theorique, note , en 

fonction des methodes de calcul choisies, et de 

1' indice deduit par interpolation de sa valeur en X±, 
i allant de 1 km, sur le spectre [A,min, Xmax] , 

- on determine 1 1 erreur Er ¥ ) , entre le spectre 
20 mesure et le spectre th6orique, 

- on minimise cette erreur en faisant varier la 
position des valeurs des indices inconnus et/ou 
l'epaisseur de la couche et/ou les valeurs des indices 
de refraction aux longueurs d'ondes initiales, et l'on 

2 5 obtient un spectre, 

- on ajoute des longueurs d ! ondes aux longueurs d'ondes 
initiales X x ... X^, les longueurs d'ondes ajoutees 
const ituant de nouveaux noeuds, 

- on repete le proc£de en choisissant un nombre m' de 
30 longueurs d'ondes initiales, m' etant sup^rieur k m, 
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et M' parametres initiaux, M' etant superieur a M, 
jusqu'a ce que la precision sur chaque spectre ainsi 
represents au mieux soit egale & une precision 
predef inie. 

5 II est question, ci-dessus, d'effectuer la 

caracterisation optique sur un intervalle de longueurs 
d'ondes X, en 1' occurrence 1' intervalle [tanin, kmax] , 
mais on ne sortirait pas du cadre de 1' invention en 
effectuant cette caracterisation optique 

10 - sur un intervalle de longueurs d'ondes 

inverses 1/A,, en 1' occurrence sur un intervalle 
[ (1/A,)min, (1/A,)max] , oil (l/A,)min est egal a 1/ (A,max) 
et (1/X)max a 1/ (Xmin) , 

ou sur un intervalle d' energies E (avec 

15 E = hv = hc/X ou h est la constante de Planck, c la 
vitesse de la lumiere dans le vide et v la frequence 
correspondant S X) , en 1' occurrence sur un intervalle 
[Emin, Emax] , ou Emin est egal a he/ (Amax) et Emax a 
he/ (A,min) , 

2 0 ou, plus generalement , sur un intervalle 

[Omin, amax] de valeurs d'une variable a fonction de X. 

L'homme du metier adaptera aisement la 
definition du procede objet de 1' invention que l'on 
vient de donner en fonction de la variable X, ainsi que 

25 les modes de realisation particuliers et les exemples 
qui en sont donnes par la suite, a l'une ou 1' autre des 
variables 1/X et E et, plus generalement, a la variable 
a dependant de X. 

Par example, si l'on effectue la 

30 caracterisation optique sur 1' intervalle [(1/A,}rnin, 
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(1/X)max] ( respect ivement sur 1'intervalle [Emin, 

Emax] ) , on remplacera 1' etape ou l'on choisit les 

longueurs d'ondes initiales par une etape ou l'on 

choisit des longueurs d'ondes inverses initiales 
5 (l/A,) i... (1/X) m (respect ivement des energies initiales 

Ei_E m ) et, plus generalement , par une etape ou l'on 

choisit des valeurs initiales (Xi-CC™ de (X. 

En outre, 1' invention s' applique pour 

caracteriser un spectre ou une partie de spectre* 
10 Selon un premier mode de mise en ceuvre 

particulier du procede objet de 1' invention, m est au 

moins egal a 2 . 

Selon un deux i erne mode, de mise en ceuvre 

particulier du procede objet de I 7 invention, m est egal 
15 . a 1 et l l on choisit des indices de refraction initiaux 

egaux. 

Selon un mode de realisation particulier de 
1' invention, le materiau est non absorbant et le nombre* 
M est egal a m, le coefficient d 1 extinction du materiau 
20 etant pris egal a 0. 

Selon un autre mode de realisation 
particulier de I 7 invention, 

- M est au moins £gal & 2 m, 

- on choisit en outre une loi d 1 interpolation pour le 
25 coefficient d 1 extinction du mat<§riau, 

. — pour chaque longueur d'onde initiale X±, 1 £ 1 £ m, 
on choisit en outre un coefficient d 1 extinction 
initial \c±, les longueurs d'ondes initiales etant en 
outre choisies de raaniere £ pouvoir determiner par 
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interpolation le coefficient d 1 extinction pour toute 
longueur d'onde de 1 1 intervalle [X min, X max] , 
— dans le processus d r optimisation, on minimise 
1 1 erreur en faisant varier en outre les valeurs des 
coefficients d 1 extinction aux longueurs d f ondes 
initiales et les longueurs d'ondes ajoutees sont en 
outre placees de fagon a representer au mieux le 
spectre du coefficient d' extinction du materiau. 

Dans ce cas, selon un mode de realisation 
particulier de 1' invention, m est egal a 1 et l'on 
choisit des indices de refraction initiaux egaux et des 
coefficients d 1 extinction initiaux egaux. 

Selon un mode de mise en oeuvre particulier 
du procede objet de 1' invention, la couche de materiau 
est mince, c'est-a-dire a une epaisseur inferieure a la 
longueur de coherence de la lumiere utilis^e pour les 
mesures, on choisit un paramdtre initial 
supplementaire, a savoir une epaisseur initiale de 
couche, et dans le processus d 1 optimisation on minimise 
1' erreur en faisant en outre varier la valeur de 
1 1 epaisseur de couche . 

S e 1 on un au t re mo de de mise en oeuvre 
particulier du proced£ objet de 1' invention, la couche 
de materiau est epaisse, c'est-a-dire n'est pas mince, 
et M est au plus egal a 2 m. 

Selon un autre mode de mise en oeuvre 
particulier, 1' epaisseur de la couche de materiau est 
connue avec une precision suf f isante et M est au plus 
egal a 2 m. 

Chaque loi d 1 interpolation peut Stre 
choisie parmi les lois d 1 interpolation lineaires, les 
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lois d' interpolation cubiques, les lois d f interpolation 
polyn&miales et les lois d 1 interpolation de type 
f one t ion « spline » . 

Dans 1' invention, la repartition des noeuds 
5 peut etre homogene . 

BREVE DESCRIPTION DES DESSINS 

La presente invention sera mieux comprise a 
la lecture de la description d'exemples de realisation 
donnes ci-apres, a titre purement indicatif et 

10 nullement limitatif, on faisant reference aux dessins 
annexes, sur lesquels : 

- la figure 1 est une vue schematique de 
dispositifs permettant de caracteriser une couche < 
conformement a 1' invention, 

15 - la figure 2 tuontre les variations de . 

l'indice de refraction en fonction de la longueur • 
d'onde, pour un materiau suivant une loi de Cauchy ..' 
(courbe I) et pour un materiau caracterise conformement 
a 1' invention (courbe II), et 

20 - la figure 3A ( respect ivement 3B) montre 

les variations de l'indice de refraction 

(respect ivement du coefficient d' extinction) en 
fonction de la longueur d'onde, pour un materiau 
suivant une loi a deux oscillateurs (courbe I) et pour 

25 un materiau caracterise conformement a 1' invention. 

EXPOSE DETAILLE DE MODES DE REALISATION PARTICULIERS 

L 1 invention propose une alternative aux 
procedes classiques, mentionnes plus haut . Elle permet 
30 d'allier la coherence d ! un modele de couche 
(correspondant & une loi d'indice continue et & des 
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epaisseurs physiques constantes) , a la generalite 
concernant la loi d 1 indice a trouver (comme dans la 
methode PAP). En outre, la resolution n'est limitee que 
par la resolution du spectre me sure . 
5 Dans le procede objet de l f invention, le 

spectre d f indice n* (X) est caracterise par : 
- un nombre reduit de "nceuds 11 , qui sont des points de 

coordonnees (A^n^ki), avec ni = n (Xi) et ki « k (Xi) , 

et 

10 - une loi d 1 interpolation entre les nceuds, qui peut 
§tre, par exemple, lineaire, cubique, de type 
« spline » ou polynomiale (de degre quelconque) . 

Cette loi d 1 interpolation permet de 
calculer, a partir des noeuds, les indices de refraction 
15 et les coefficients d 1 extinction pour les longueurs 
d'ondes situees entre les noeuds. 

Par exemple, lorsque 1 1 indice de refraction 
est caracterise par un ensemble de valeurs pour des 
longueurs d'ondes X x ... X^, on peut utiliser une 
20 interpolation lineaire entre deux longueurs d ! ondes X^ 
et Xi+i pour calculer 1 1 indice n a la longueur d'onde X 
(voir le document [6] chapitre 3) : 

n (X) - -*M4) + (^-4)"(4+i) 
avec A,i<X<A,i + i 

25 On peut f aire de meme pour le coefficient 

d' extinction. 

Lorsque le nombre de nceuds le permet, des 
formules d 1 interpolation plus complexes, faisant 
intervenir les nuauds voisins, peuvent etre utilisees 

3 0 (voir le document [6] chapitre 3) . 
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Un modele de couche est done caracteris^ 
par une gpaisseur 8 et une famille de noeuds. 

On donne ci-apres un exemple du procede 
objet de l 1 invention. Dans cet exemple, les mesures *F 
5 sont constitutes d'une mesure ref lectometrigue R(X) et 
d'une mesure ellipsometrique Si, 2 (9A) oH Q est 1' angle 
d' incidence du faisceau lumineux que l'on envoie sur la 
couche a etudier lors de la mesure d f ellipsometrie . 

Cette couche est une couche mince de sorte 
10 que l 1 epaisseur de cette couche est aussi une variable 
du probleme. En outre, on suppose qu ! une seule couche 
est inconnue , cette couche etant f ormee sur un substrat , 
connu . 

Expliquons d f abord bridvement cet exemple, 
15 . qui utilise un algorithme. 

A partir d 1 informations supposees sur 
l 1 epaisseur £ de la couche etudiee, sur 1 1 indice de 
refraction n(X) et sur le coefficient d ? extinction 
k(X), du mat£riau de cette couche, on construit des 
2 0 noeuds de depart (en faible nombre) et une epaisseur de 
depart s. 

On dispose ainsi de m noeuds et, par 
interpolation, on peut connaitre n{\) et k(X) en dehors 
des valeurs des longueurs d'ondes associees aux noeuds. 
25 A partir de 1' epaisseur de depart £ et de 

ces valeurs de depart n(X) et k(X), on determine le 
spectre theorique en utilisant des calculs 

ellipsometriques et ref lectrometriques . 

Par ailleurs, au moyen de dispositifs 
30 d' ellipsometrie et de ref lectrometrie et d'un 
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spectrometre, on obtient Si, 2 (9 /A,) et R{\) et 1 1 on en 
deduit les mesures notees *F . 

On compare ensuite ¥ et ¥ en utilisant une 
fonction d'erreur Er et l'on optimise la valeur de 
5 1 r indice de refraction et la valeur du coefficient 
d 1 extinction aux diff brents noeuds, ainsi que la valeur 
de l 1 epaisseur, en cherchant a minimiser Er ( % ) . 

Lorsque ces valeurs sont optimisees et si 
la precision sur le spectre n(X) , le spectre k(k) et 
10 1 1 epaisseur £ n 1 est pas suf f isante, on ajoute de 
nouveaux noeuds, on fait varier I 1 epaisseur £, et l'on 
recommence la determination de la comparaison de *P 

et *F et 1 1 optimisation raentionnees plus haut, etc. 

On arr§te la boucle ainsi definie lorsque 

15 la precision sur chacun des spectres n(X) et k(A,) et 

sur 1' epaisseur 8 est jugee suf f isante. 

Les spectres n(X), k(X) et I 1 epaisseur e 
sont ainsi caracterises . 

La figure 1 montre de fagon schematique la 

2 0 couche etudiee 2, formee sur un subs t rat 4. On voit le 

dispositif d r ellipsometrie 5,6, le dispositif de 
reflect rometrie 8 et le spectrometre 10. On voit en 
outre des moyens electroniques de traitement 12 
permettant de caracteriser n(A.) , k(X) et £ en fonction 
25 des informations fournies par le spectrometre 10 et 
conformement au procede de 1' invention. Ces moyens 12 
sont munis de moyens d'affichage 14. 

Revenons de faqron plus detaillee sur 
1 1 exemple donne , 

3 0 Phase 1 
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Le procede de cet exemple comporte d'abord 
une etape d 1 initialisation. 

L 1 algorithme debute avec tin nombre reduit 
de noeuds (2 ou plus de 2) . Cependant, on peut debuter 
par un seul noeud en imposant un indice de refraction et 
un coefficient d ! extinction qui restent constants 
lorsque la longueur d ! onde varie . 

On choisit des positions de noeuds de fagon 
a pouvoir, a partir de cette famille de noeuds, deduire 
tout le spectre par interpolation. Le modele de couche 
est done a 3 parametres ou plus, puisque l'epaisseur 
est aussi une variable a determiner. La table d 1 indice 
sur tout le spectre est done d^duite des noeuds par 
interpolation . 

C ! est le cas lorsque l'epaisseur des 
couches est de 1 1 ordre de grandeur de la longueur 
d f onde (couches minces). Mais lorsque l'epaisseur de la 
couche est tres importante, l'epaisseur n'intervient 
quasiment plus dans le calcul de la reponse de la 
couche, et n'est done plus une variable du probleme . 

Par exemple, dans le cas de disques 
optiques, sur lesquels une couche tres import ante est 
deposee (son <§paisseur est de 1 ' ordre de 1 millimetre), 
le coefficient de reflexion d'une telle couche n'est 
plus fonction de l'epaisseur de la couche, mais 
seulement de 1 f indice de refraction de celle-ci. 

Par exemple, on choisit de placer les deux 
premiers noeuds aux extremites X min et X max du 
spectre. Les valeurs de 1 1 indice complexe en ces 
extremites sont choisies en fonction du type du 
mat6riau etudie. Par exemple, sur un spectre 
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ellipsometrique entre 3 00nm et 800nm d'une couche mince 
de resine photosensible, on prenci n (3 00nm) =n (800nm) =1,5 
et k(300nm) =k(800nm) =0 . 

Lorsque le spectre n'est caracterise que 
5 par deux ntEuds, 1 1 indice entre les extremes est 
determine par interpolation lineaire. Dans le cas 
considere, on a done n(A,)=l,5 et k(A,)=0 pour X 
appar tenant a [3 00nm, 800nm] . 

A partir de trois noeuds, on choisit plut6t 
10 une interpolation cubique afin d'obtenir des formes de 
loi d f indice plus douces que celles obtenues par 
interpolation linSaire. 

L'epaisseur de depart est, quant a elle, 
choisie aussi proche que possible de I'gpaisseur 
15 reelle. 

Phase 2 

On procede ensuite a une determination 
optimale des valeurs de l 1 indice de refraction et du 
coefficient d 1 extinction sur les nosuds et de la valeur 
20 de l'epaisseur. 

Pour ce faire, les spectres *F (X) sont 
calcules a 1 1 aide du moddle de couche utilise, 
resultant du choix des nceuds, de la loi d 1 interpolation 
et de l'epaisseur de la couche. 

25 Le modele physique utilise pour le calcul 

de X F est bien sur f one t ion de la methode de me sure 
utilisee , c 1 est-a-dire notamment de 1 1 angle d f incidence 
de la liimiere, du spectre utilise et du modele de 
couches minces ou de couches epaisses le cas echeant 

30 (voir par example les modeles de couches empilees dans 
le document [3 J ) . 
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Les spectres *F etant constitues par un 
ensemble de mesures de natures diverses (par exemple 
des mesures ellipsometriques et ref lectometriques) , on 
utilise pour les mesures de ref lectometrie 
( respect ivement d 1 ellipsometrie) uhe methode de calcul 
de rSf lectometrie (respect ivement d 1 ellipsometrie) . 

Les mesures ref lectometriques et 

ellipsometriques sont combinees par 1 1 intermediaire 
d ! une fonction d'erreur Er( x i / / T F) qui est par exemple du 
genre de celle qui est definie par 1* equation (1) - 

GrSce a line fonction de recherche telle que 
Simplex par exemple, on minimise l'ecart entre *F {X) et 
"¥(X) 9 en faisant varier la valeur de 1 1 indice de 
refraction et la valeur du coefficient d f extinction a 
la position de chacun des nceuds ainsi que l 1 epaisseur 
de couche (si cette epaisseur est un facteur important 
dans le calcul de *F ) . Lorsque 1 1 ecart est minimum 
c'est-a-dire lorsque Er (¥,*?) est minimum, cela 
signifie que les mesures ref lectometriques et 
ellipsometriques coincident au mieux (pour un nombre de 
nceuds donne ) . 

A ce stade, on obtient, pour un nombre 
connu de noeuds et une position spectrale connue pour 
chacun de ces des nceuds, le modele de couche (indice de 
refraction, coefficient d 1 extinction et epaisseur) qui 
correspond au mieux a la couche reelle. 

La validity du modele trouv6 est d'autant 
plus assuree que le nombre de mesures est grand. Pour 
avoir un grand nombre de mesures, on peut par exemple 
utiliser plusieurs angles d 1 incidence de la lumiere &±, 
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1 £ ± £ £ , pour les spectres ellipsometriques, faire 
une mesure ref lectometrique et faire des mesures 
goniometriques complement aires . 
Phase 3 

5 Ensuite, on ajoute un nombre fini de noeuds . 

Dans un premier mode de realisation, les noeuds rajoutes 
sont positionnes de fagon a representer au mieux les 
spectres n(X) et k(A,) . A titre d'exemple, on place ces 
noeuds supplementaires aux endroits ou 1 1 ecart entre *F 

10 et est maximum ou aux endroits en lesquels les noeuds 
sont les plus espaces . Et l ! on retourne la phase 2 
tant que la precision sur chacun des spectres xi(X) et 
k(A,) et sur I'epaisseur 8 n'est pas suffisante, 
c r est-a-dire n'est pas egale a une precision 

15 predefinie. 

Dans un deuxieme mode de realisation, on 
intercalle de nouveaux noeuds entre deux noeuds de 
1' ensemble de noeuds prealablement choisis et qui se 
r^partiront de fagon homog^ne sur le spectre. 

2 0 Dans ce qui suit, nous donnons deux 

exemples courants d 1 application de 1 1 invention . Ces 
exemples met tent en jeu deux types de materiaux qui 
suivent des lois differentes. 

A partir de mesures ellipsometriques et 
25 ref lectometriques, nous proposons de retrouver les lois 
physiques suivies par ces materiaux. Nous procedons de 
la fac^on suivante : 

Un materiau fictif est cree, ce materiau 
suivant une loi theorique connue (une loi de Cauchy ou 

3 0 une loi d 1 oscillateurs harmoniques) , avec des 

parametres que nous fixons arbitrairement . Les 
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variations de 1 1 indice de refraction complexe en 
fonction de la longueur d'onde sont ainsi parf aitement 
connues . 

Nous imposons en outre au materiau une 
epaisseur de 200,00nm, sur un substrat de silicium, ce 
dernier etant trds bien connu. 

Des mesures fictives (mesures 

ellipsometriques, mesures r€f lectometriques) sont 
calculees, puis bruitees de fa<?on a introduire un 
d#faut d* appareillage. 

Tout se passe comme si nous avions des 
mesures r<§elles, effectuees sur le materiau. Mais, 
contrairement a la realite, nous connaissons 
parf aitement l 1 indice de refraction complexe puisque 
nous l'avons fixe, de meme que nous avons fixe 
1* epaisseur de la couche de matSriau. 

Nous testons ici la methode "en aveugle", 
c'est-&-dire que nous partons d'une epaisseur fausse 
(220nm) et d ! indices de refraction complexes faux, 
puisqu ! ils sont senses §tre inconnus . 

Nous appliquons le proc^de objet de 
1 1 invention puis nous comparons 1 1 indice de refraction 
complexe trouve a 1 1 indice de refraction complexe 
theorique . Nous retrouvons bien les memes lois , tres 
precisement, ainsi que la meme epaisseur de couche. 

Prenons comme premier exemple un materiau 
dont 1' indice de refraction complexe suit une loi de 
Cauchy telle que : 

ir ai300 2 A1 300 4 
n(X) = 1,5 + 0,1-^- + <U-^r 

k(X)=0 
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Cette loi d'indice est typique des resines 
photosensibles (dans la plage spectrale 300-800nm) . 

Afin de retrouver a 1 1 aide de la methode 
des noeuds (c'est-a-dire du procede objet de 
5 1 1 invention) la loi d ! indice mentionnee ci-dessus, nous 
effectuons deux mesures, a savoir vine tnesure 
ellipsometrique a un angle de 70° et une mesure 
ref lectometrique . 

Les conditions de traitement sont les 

10 suivantes : 

- le spectre traits est compris entre 300nm et 800nm, 

- au depart, les noeuds sont aux positions (300nm, 1,6) 
et (800nm, 1,6), c'est-a-dire que l'indice est 
considere comme variant lineairement entre 3 00nm et 

15 800nm, et que sa valeur est constante (6gale a 1/6), 

- le notnbre de noeuds est augmente it€rativement en 
rajoutant au moins un noeud a chaque iteration (par 
exemple, on ajoute 2 nceuds) , 

- la loi d f interpolation est une loi cubique, lorsque 
20 le nombre de noeuds est superieur a 2, sinon elle est 

lineaire, et 

- 1 1 algorithme de minimisation utilise est un 
algorithme de type Simplex. 

La figure 2 permet de comparer l'indice de 
25 refraction correspondant au materiau fictif qui suit 
parfaitement une loi de Cauchy (courbe I) a l^ndice de 
refraction que nous trouvons par la methode des nceuds 
(courbe II), a I 1 aide de 6 noeuds (representes par des 
cercles sur la figure 2) . Nous avons utilise un 
30 compose d'une mesure ellipsometrique {S x (X) , S 2 (X) } a 70° 
et d'une mesure ref lectometrique R(X) . 
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Prenons comme deuxieme exemple un materiau 
dont l'indice de refraction complexe suit une loi a 
deux oscillateurs, telle que : 

\n(E) + jk(Ejf = 1 + V ^ ^ 

5 avec 

j 2 = -1 et E = 1240/A, en nm) 
Ai=0 , 25x1 , 5 2 A 2 =0 , 25x1 , 5 2 
Ei=124 0/400 E 2 =124 0/300 
Gi=0,3 G 2 =0,3 
10 La methode des noeuds est appliquee a un 

ensemble de mesures ellipsom^triques effectuees entre 
250nm et 800nm, h 75° , 70° , 60° , et 45°, avec en plus 
une mesure ref lectometrique . L'epaisseur r§elle du 
materiau etant de 200nm / on trouve une epaisseur de 
15 2 00,25 nm avec la methode des nceuds . L'ajustement sur 
la loi d'indice est trls bon. 

Les figures 3 A et 3B montrent 
respectivement les reconstructions des courbes n(X) et 
k(A,) du materiau a l T aide de la methode des noeuds. La 
20 reconstruction est effectuee a I'aide de quatre 
spectres ellipsometriques et d'un spectre 

r£f lectometrique . Les pics d 1 absorption reels sont tres 

bien repr§sentes par la courbe obtenue par 
* 

interpolation cubique entre les noeuds (represents par 
25 des cercles sur les figures 3A et'3B) . 

Sur la figure 3A, la courbe I 
(respectivement II) correspond a n suivant parfaitement 
la loi choisie (respectivement & n trouve par la 
methode des noeuds) . 
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Sur la figure 3B, la courbe I 
(respectivement II) correspond a k suivant parf aitement 
la loi choisie (respectivement a k trouve par la 
m^thode des nceuds) . 
5 On v i ent de de c r i r e de s exerap 1 e s de 

1' invention. On remarquera d'une maniere plus generale 
que, dans cette derniere, on considere un ensemble X de 
valeur s , avec X= {n x , n 2 , ... , n± , ...n m , ki , k 2 , ... , ki / ... / k m , 8} , ou 

m est la valeur de 1 1 indice de refraction 
10 au nceud X±, ie{l...m}, m etant le norabre de noeuds 

k± est la valeur du coefficient d 1 absorption 

au noeud X±, ie{l...m}, 

E est l'epaisseur de la couche etudi<§e. 
L 1 operation de minimisation d'erreur Er (¥, 
15 W) revient a trouver 1' ensemble X tel que Er soit 
minimum . 

Lorsque 1 1 on n' impose pas de contrainte 
particuliere, la minimisation est une minimisation a 
2xm+l parametres. On peut bien sur a j outer des 
20 contraintes afin de diminuer le nombre de variables. 

En particulier, si 1' on sait que le 
materiau est non absorbant, on impose ki=0 pour tout i 
de {l...m} et X devient : X={n x , n 2 ..., ni , ...n m/ £} . 

Si, par une mesure complementaire (par 
25 exemple tine mesure de goniometrie ou une mesure directe 
non optique) , on connait l'epaisseur de la couche 
consideree avec une precision suffisante, l'epaisseur 8 
n'est plus une variable et I 1 on a : 

X^fni , n 2 , ... , Hi , ...n m , ki , k 2 , ... , ki , ... , k m } . 
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Les deux options precedent es peuvent £tre 
bien sur corabinees . 

De plus, la presente invention n'est pas 
limitee Sl la caracterisation de couches minces. Elle 
s' applique aussi & la caracterisation de couches 
€paisses. 

En outre, la presente invention n'est pas 
limitee a la caracterisation d'une seule couche, forrmie 
sur un substrata Elle s ' applique aussi a la 
caracterisation de deux, ou plus de deux, couches 
formSes sur un substrat . 

Les documents cit£s dans la presente 
description sont les suivants : 

[1] R.M.A. Azzam and N.M. Bashara, 
Ellipsometry and Polarized Light, North-Holland Physics 
Publishing, 1997, chapitre 6. 

[2] B. K. Minhas, S.A. Coulombe, S. Sohail, 
H . Naqvi and J.R. McNeil, Ellipsometric scatterometry 
for the metrology of sub-0 . l-p,m-linewidth structures,' 
Applied Optics, 37(22) : 5112-5115, 1998. 

[3] M. Borri and E . Wolf, Principle of 
Optics, Cambridge University Press edition. 

[4] A.R. Forouhi and I. Bloomer, Optical 
dispersion relations for amorphous semiconductors and 
amorphous dielectrics, Physical Review B, 34(10) : 
7018-7026, November 1986. 

[5] F. L. Terry, Jr., A modified harmonic 
oscillator approximation scheme for the dielectric 
constants of Al x Gi- x As, Journal of Applied Physics, 
70(1), 1991, pages 409-417. 
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[6] W. H. Press, S. A. Teukolsky, W. T. 
Vetterling and B. P. Flannery, Numerical Recipes in C, 
Cambridge University Press, 1992, chapitres 3 et 10. 
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RE VEND I CAT I ONS 

1. Procgde de caract6risation optique d'au 
moins une couche d'un materiau dans un intervalle de 
longueurs d'ondes [X min, X max], cette couche etant 
5 form§e sur un substrat, ce proc^de 6tant caracterise en 
ce que : 

- on effectue un ensemble de mesures de ref lectomStrie 
et/ou d r ellipsometrie, cet ensemble de mesures 
conduisant a un spectre mesure, note *F, 

10 - on choisit m longueurs d'ondes initiales X x •«X ro 
appartenant a cet intervalle, m etant un nombre entier 
au moins egal & 1, on associe, a chaque longueur 
d'onde, un indice de refraction, 

- on choisit une loi d 1 interpolation au moins pour 
15 I 1 indice de refraction du materiau, pour les longueurs 

d'ondes comprises entre les longueurs d'ondes 
initiales Xi-Am, 

- on choisit M paramdtres initiaux, M etant au moins 
egal §. m, a savoir un indice de refraction initial ni 

20 pour chaque longueur d'onde initiale Xj., 1 ^ i ^ m, 
les longueurs d'ondes initiales etant choisies de 
man i ere a pouvoir determiner par interpolation au 
moins I 1 indice de refraction pour toute longueur 
d'onde de 1 ' intervalle [X min, X max], les couples 

25 (A,j.,ni) etant appeles nceuds, 

- on choisit des methodes de calcul de r<§f lectometrie 
et d ! ellipsometrie, 

- on choisit aussi une fonction d'erreur Er, 
representative de l'ecart entre deux spectres et 

3 0 les spectres *¥ x et *F a 6tant calcules ou mesures 
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sur un nombre de points super ieur au nombre m de 
n<Euds , 

— a l ! aide des m longueurs d f ondes initiales, des M 
parametres initiaux et de la loi d* interpolation, on 

5 met en aeuvre le processus d 1 optimisation suivant : 

— on determine un spectre theorique, note 

*P , en fonction des methodes de calcul choisies, et de 
1'indice deduit par interpolation de sa valeur en Xi, i 
allant de 1 a m, sur le spectre [A,min, A,max] , 

10 — on determine 1 1 erreur Er (*P , ) , entre le spectre 
mesure et le spectre theorique 

— on minimise cette erreur en faisant varier la 
position des valeurs des indices inconnus et/ou 
l'epaisseur de la couche et/ou les valeurs des indices 

15 de refraction aux longueurs d'ondes initiales, et 1'on 
obtient un spectre, 

— on ajoute des longueurs d'ondes aux longueurs d ! ondes 
initiales Ai X^, les longueurs d'ondes ajoutees 
constituant de nouveaux nceuds, 

20 — on repete le procede en choisissant un nombre m' de 
longueurs d'ondes initiales, m' efcant superieur a m, 
et M' parametres initiaux, M' etant superieur a M, 
jusqu f a ce que la precision sur chaque spectre ainsi 
represents au mieux soit egale a une precision 

25 predef inie . 

2. Procede selon la revendication 1, dans 
lequel m est au moins egal a 2 . 

3. Procede selon la revendication 1, dans 
lequel m est egal a 1 et 1 1 on choisit des indices de 

3 0 refraction initiaux egaux. 
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4. Proced6 selon l'une quelconque des 
revendications 1 a 3 , dans lequel le materiau est non 
absorbant et le nombre M est egal a m, le coefficient 
d' extinction du materiau etant pris £gal a 0. 
5 5. Procede selon l'une quelconque des 

revendications 1 a 3, dans lequel : 

— M est au raoins egal a 2 m, 

— on choisit en outre une loi d' interpolation pour le 
coefficient d r extinction du materiau, 

10 - pour chaque longueur d'onde initiale X±, 1 £ i £ m, 
on choisit en outre un coefficient d' extinction 
initial ki, les longueurs d'ondes initiales etant en 
outre choisies de manilre & pouvoir determiner par 
interpolation le coefficient d 1 extinction pour toute 

15 longueur d'onde de 1 1 intervalle [X min, X max] , 

— dans le processus d 1 optimisation, on minimise 
1 ' erreur en faisant varier en outre les valeurs des 
coefficients d' extinction aux longueurs d'ondes 
initiales et les longueurs d f ondes ajoutees sorit en 

20 outre placees de fagon a representer au mieux le 

spectre du coefficient d 1 extinction du materiau. 

6. Procede selon la revendication 5, dans 
lequel m est egal a 1 et 1 ' on choisit des indices de 
refraction initiaux egaux et des coefficients 

25 d ! extinction initiaux egaux . 

7. Procede selon l'une quelconque des 
revendications 1^6, dans lequel la couche de materiau 
est mince, c'est-&-dire a une epaisseur inf6rieure a la 
longueur de coherence de la lumiere utilis^e pour les 

3 0 mesures, on choisit un parametre initial 
supplement aire, a savoir une Epaisseur initiale de 
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couche, et dans le processus d 1 optimisation on minimise 
l'erreur en faisant en outre varier la valeur de 
l'epaisseur de couche . 

8. Procede selon I'une quelconque des 
5 revendications 1 a 6, dans leguel la couche de materiau 

est epaisse, c'est-a-dire n'est pas mince, et M est au 
plus egal a 2 m. 

9. Procede selon 1 1 une quelconque des 
revendications 1 a 6, dans lequel l'epaisseur de la 

10 couche de materiau est connue avec une precision 
predefinie et M est au plus egal a 2 m. 

10. Proc£d£ selon I'une quelconque des 
revendications 1 a 9, dans lequel chaque loi 
d' interpolation est choisie parmi les lois 

15 d ! interpolation lineaires, les lois d 1 interpolation 
cubiques, les lois d 1 interpolation polynomiales et les 
lois d 1 interpolation de type fonction spline . 

11. Procede selon I'une quelconque des 
revendications 1 a 10 , dans lequel la repartition des 

20 nosuds est homo gene . 
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